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^ Abstract. We present a perturbative result for the temporal evolution of the fidelity 

of the quantum kicked rotor, i.e. the overlap of the same initial state evolved with two 
slightly different kicking strengths, for kicking periods close to a principal quantum 
resonance. Based on a pendulum approximation we describe the fidelity for rotational 
orbits in the pseudo-classical phase space of a corresponding classical map. Our results 
, are compared to numerical simulations indicating the range of applicability of our 

analytical approximation. 
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(N ! 1. Introduction 

> 

O ■ In classical mechanics chaos can be defined using the stability of a trajectory. Consider 



two neighbouring phase space points, the distance of their orbits in the phase space will 
grow exponentially in time for a generic chaotic system. Such an approach is bound 
O ■ to fail in a quantum mechanical treatment as the time evolution is unitary and the 
overlap of two wave packets is constant in time. However changing a parameter of the 
Hamiltonian instead of the initial state will lead to an overlap varying in time. This idea 
■ was formulated by Peres [1] when he introduced the fidelity, also known as Loschmidt 
^ ■ echo |2j- The fidelity is defined as the overlap of an initial state evolved with slightly 
different Hamiltonians. This quantity has been used to characterise the stability of 
quantum states [3J. For classically chaotic systems, the fidelity shows some generic 
behaviour [3]. 

Here we study the fidelity in the quasi-integrable regime for which not so many 
results are established jlHH]. We will do this for one of the best known examples of 
classically chaotic system: the kicked rotor (KR) [9l[T0]. In order to understand the 
fingerprints of classical chaos in quantum mechanics the quantum kicked rotor (QKR) 
has been and still is a fruitful field of study pjj. It shows several interesting phenomena 
like quantum resonance [12] and dynamical localisation |13] both in direct contradiction 
to the behaviour of the classical system. In the implementation of the QKR near a 
quantum resonance in the gravity field Oberthaler et al. observed quantum accelerator 
modes [H]. Fishman et al. were able to describe these theoretically using a pseudo- 
classical limit. There the detuning of the kicking period to its resonant value plays the 
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role of the Planck constant [TBj. Using this treatment the QKR can be mapped onto a 
kicked rotor with a renormalized kicking strength. This leads to regular structures and 
allows the application of semi-classical methods in the pseudo-classical limit although 
the system might be chaotic in the true semi-classical limit |T6] . 

The KR shows two types of motion in the quasi-integrable regime, oscillations about 
the stable fixed points and rotating motion. The motion on the classical resonance island 
surrounding the fixed point leads to revivals of the quantum fidelity. This was shown for 
the regular QKR by Sankaranarayanan et al. |5] and for the near resonant QKR by Abb 
et al. j6]. Similar results were obtained for another kicked system by Krivolapov et al. [7]. 
Rotational and oscillating orbits were also numerically studied in |[8j. Our focus in this 
paper is to treat the rotating modes of the QKR near a quantum resonance. Therefore 
we will use the pseudo-classical method and apply the pendulum approximation. Section 
[2] sets the stage by reviewing the pseudo-classical approximation and defining the fidelity 
which is studied here. In Section|3]we will give a perturbative treatment of the pendulum 
and discuss its validity based on a comparison with numerical simulations in Section HI 
An additional numerical check is shown in Sect. 5. 

2. Fidelity for the atom optical kicked rotor 

The story of the experimental investigation of the quantum kicked rotor is quite long [17] 
and it has been continuing until today (see, e.g., the refs. [T81 - I20] and, for fidelity 
measurements specifically, the refs. [21]). In the experimental realisations, atoms 
are kicked by a periodic potential formed by a standing wave of laser light, i.e. an 
optical lattice which is flashed on and off periodically in time [T71|22]. Using rescaled 
dimensionless momentum p, position x, kicking period r and kicking strength k, see 
e.g. [15], the Hamiltonian for one kicked atom is: 

n OO 

H{p.,x,t) = — + kcos{x) 6{t — nT). (1) 

n=— oo 

The Floquet operator mapping the state right after one kick to the state right after the 
next kick is, see e.g. [TTj . 

^ ^ ^-ikcos{X)^-i^P\ (2) 

The dynamics is obtained by repeated application of this operator. In contrast to the 
usual kicked rotor, a kicked atom lives along a line. Doing a gauge transformation one 
can still recover a problem with conserved quasi-momentum, see e.g. [T5|[T6]. Here the 
quasi- momentum corresponds to the fractional part of the momentum p = n + P, where 
n is an integer and /3 is a real variable between and 1. The problem of one atom 
along a line is mapped onto a problem of a continuous family of rotors. Each of them 
corresponds to one value of /3. That is the reason why we will speak from now on of a 
/3— rotor. In realisations of the kicked rotor with Bose-Einstein condensates it has been 
checked that the interactions between atoms in the cloud can be neglected [IHlllS] . This 
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brings another justification of our one particle approacli. Tlie wave function of a single 
/3— rotor |\E'^) is obtained from the wave function of the kicked atom by p^[T6]: 



while the Floquet operator for one /9— rotor is 

^^^g-ifccose~g-i^(A^+/3)2^ (4) 

where J\f is the angular momentum operator. The operator (jlj) formally differs from the 
usual Floquet operator [llj only by the quasi- momentum f3. 

A principal quantum resonance occurs in the kicked particle when the phase due 
to the free evolution in vanishes [23]. Each resonance leads to a ballistic motion 
and a quadratic growth of the energy. This happens at r = 27rZ and for resonant quasi- 
momenta Pres = 1/2 + q/ 1 for integers / and q such that / > 1 and < g < / — 1. To 
simplify the notation, in this paper we consider the specific resonance with t = 2tt and 

/Sres = 1/2. 

Take now a kicking period slightly detuned from its resonant value r = 27r + e. 
Introducing the rescaled momentum / = \e\J\f we can rewrite the Floquet operator (j4]) 



as 



This operator can be identified as the formal quantisation of another fictitious kicked 
rotor. The main benefit of this mapping is that the kicking strength of the new problem 
is k = \e\k and |e| plays the role of the Planck constant. We will be interested from now 
on in the regime |e| — )■ 0, called e— semiclassical limit [15J, which must not be confused 
with the true semiclassical one of the initial problem. In the e— semiclassical limit one 
can derive easily the e— classical map, which is very similar to the celebrated standard 
map [9]: 

h+i = h + ksiYi{Ot+i) 

e^^^ = et + sgn(e)/t -71 + T/3 mod 27r . (6) 

We are interested in the dynamics around a stable fixed point of One common 
approximation is the pendulum approximation [9|[T0]: 

Hp,M,d)=^—^ + kcose, (7) 

where we have defined ^ = sgn(e)(— vr + r/3). For a given (even large) kicking strength 
k, one can always choose a small enough e so that we are in the quasi-integrable regime. 
This makes the phase spaces corresponding to (EI) and (^^, respectively, very similar, see 

Fig.m 

The main goal of this paper is to check the stability of the quantum dynamics 
under slight variation of the kicking strength. The fidelity |Tl|3] appears as the natural 
quantity to look at. For one single /3— rotor it is defined as: 

Fpik„k2,t)= i^PoKl^i^L^m'. (8) 
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Figure 1. (Colour online) Phasespace of the pendulum (black ) and the kicked 

rotor (grey ) for k — \e\k = O.OStt. 



For the initial problem of a kicked atom one needs to consider the fidelity for a sub- 
ensemble of rotors, which is defined then by |M1[25] : 

" -A/3 2 



F{ki,k2,/3i,A/3,t) 



d/3 {^Ull^UlJ^o) 



(9) 



If the initial state mainly lives on the stable island the fidelity shows revivals as explained 
in [6], see also a similar context in [7]. Here the study will be devoted to the fidelity ([8]) 
in the neighbourhood of such an island. We will approximate (IH]) by the fidelity of the 
pendulum (^^. The fidelity of the pendulum is obtained by expanding the initial state 
in the eigenbasis of the Hamiltoniaiiil of Eq. ([7]), which depends on /3 via ^. 



Fp{ki,k2,t) = 

Y.i'^^it = 0)\Mk2)){Mk2)\Mki)){Mki)\'^pit = 0))e' 



(10) 



Throughout the paper we focus on the quantum problem associated to ([5]) so that |e| 
is always our Planck constant. This is the reason why we will call the e— semiclassical 
regime simply the semiclassical regime. 



3. Perturbative treatment of the pendulum 

We are interested in the quantum pendulum following the approximation ([7]). It is well 
known that this system is quantum mechanically integrable [26]. Following ( ITOl) we 
want to obtain simple explicit formulae for the eigenenergies and the eigenfunctions 
\<pn{k)) of the Hamiltonian when |e| is going to 0. We will follow standard perturbation 

X We emphasise the dependence of the eigenstates on k as it is the perturbation parameter in the 
fidelity. 
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theory. The only unusual thing is that the potential is proportional to the effective 
Planck constant so that it vanishes at the classical limit |e| = 0. 

Our unperturbed system is a free particle along a ring with eigenenergy and 
eigenfunction (m G Z): 

En. =^^^^^^ = |+emH + 0(|er), (11) 

„im8 



For A; > the Schrodinger equation for the stationary states becomes: 

-i|e|— + n ^ + kcose^! = E^! . (13) 



2 V 9^ 

Doing the gauge transformation \& = exp(— i^6'/|e|)^/' and setting ip{6) = f{z = 9/2), 
Eq. ( fT3l) becomes: 

dV f 8E Ak , ,\ ^, , 

^+(^-2^cos(2^)j/(z) = 0, (14) 

which is the standard form of Mathieu equation, see e.g. 16.2.1 p. 97 in [27]. For our 
purpose it is easier to look for a solution of f|T^ as the following series: 

/(z) = ^c.e(^+2in). ^15) 

As we require the "true" wave function "^{6) to be univalued one can easily see that we 
need 

/^=^- (16) 



The eigenenergies of f[T3|) are given by characteristic values of Mathieu functions, which 
do not lead to simple explicit formulae. A semiclassical approach is rather used to write 
an expansion of the eigenenergies. The details are found in the Appendix A The results 
are, assuming k = k\e\ and noting = sgn(e)7r(2/3 — 1): 



3 (m + ^)2fc2 ^ 5 k^ 

?4 



50 



4 a ^sisji'i- 

In (fT5|) the coefficients c„ are the solutions of the following recurrence relation: 

[2E - {n\e\ + 0'] c„ = ~k{cn-i + c„+i) . (18) 

If we assume now that we start from an unperturbed state (IT^ with the energy (ITT]) . 
([T8|l can be rewritten as: 

(n-m)eo (^) _ (^) (^) 
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which gives the solution: ci™'' = Jn~m{—k/C,o) where J„(x) stands for the Bessel function 
of integer order n. The perturbed eigenf unctions are then: 



{9\(PUk)) = y,Jn-m (^]e 



in6 gimS— ifcsin(6)/5o (20) 



where we have used in the second equahty the following identity for the Bessel functions, 
see e.g. 7.2.4(26) p.7 in [28]: 

^ J„(x)e'"^ = e'"^'°V (21) 

The great benefit from (120|) is that we can directly compute the overlap coefficient for 
the fidelity ( ITOj) . We assume that the initial state is a plane wave with momentum uq. 
{nl'^isit = 0)) = 6n,no where 6n,k is the Kronecker symbol. Then: 

(0^(A;)|vl/^(t = 0)) = J2{'Pmmn){n\Mt = 0)) = J„o— (^) , (22) 

(0m(fcl)|0n(fc2)) = J^(0m(A;i) b) (p|0n(fc2)) 



Jm—n I >. 1 • (23) 



^0 

Finally our simple perturbative approach lets us write an explicit formula for the fidelity, 
reminding = sgn(e)7r(2/3 — 1): 

Fp{k^,k2,t) = 

^ ^ Jm-n ( Z I Jno-m ( ~~Z ) Jno~n ( ~7 ) ^ , (24) 

' n.mez V 4o y V 4o y V 4o y 

where i?^^ and E^^ are given by f|T7|) . The formula f l24|) is the main result of this paper. 
In the next section we will estimate numerically the range of validity and the accuracy 
of to describe the quantum kicked rotor. In Fig. |2]we can already see that the more 
orders we take for the energy, the more accurate we get. 



4. Numerical comparison of the approaches 

The perturbative approach, cf. Eq. f l2^ . will now be checked numerically in this section. 
As single rotors and ensembles show qualitatively different behaviour we treat these cases 
separately. 

Before showing the main results, we discuss the intrinsic limitations of our approach. 
The phase space of the pendulum is the cylinder whereas the phase space of the KR 
can be mapped onto a torus. The pendulum approximation is therefore only valid in 
one phase space cell of the KR. Due to this mismatch we expect the approximation to 
fail at the border of a cell. 

In the derivation of the perturbative result we explicitly focused on the rotating 
regime, which means that our results have to fail in the description of states on the 




Figure 2. (Colour online) Fidelity using the pendulum (green ), the perturbative 

result with the third (black ) and fourth (blue — • • — ) order in |e| in the energy, 

and the original QKR (red ). /? = 0.3, e = 0.05, ki = O.Gtt and = O.Stt. The 

data are averaged over 100 kicks in order to cancel fast oscillations. 



island. The half width of the island in the pendulum approximation is given by 
AI = 2^/k\7\ [in], which becomes in units of the quasi-momentum 

AP«. = ^.M. (25) 

T V TT"^ 

When the distance from /3 to its resonant value (centred at the island) is less than A/3c 
we expect to observe the revivals in the fidelity as described in [6]. 

4-1- Single rotors 

The most important goal of the discussion of single rotors is to get an intuition for the 
quality of the pendulum approximation. We can discuss this step only for single rotors 
as the calculation of the pendulum fidelity is numerically very challenging. Averaging 
the fidelity over 100 kicks allows to identify maxima in the two cases and to read of the 
amplitude and the period (c.f., for instance, Fig. |3]). Plotting the relative deviation of 
the period of the maxima one observes that this relative error is nearly independent of 
the choice of the maximum. As a measure of deviation of the amplitude we compared 
a limited number of maxima, whilst these maxima should be visible both in the QKR 
and the pendulum data. We decide to take the maximal deviation within the first 10 
maxima. 

When requiring an accuracy of 10% we can give as a boundary of validity f3 > 0.2. 
In the case of the amplitude the criterion is not as clear. Taking also the deviation in 
the amplitude into account we concluded that (3 > 0.3 following this criterion. 

4-2. Ensembles 

In order to build ensembles we need to evaluate the integral in (Q. This integral 
is approximated by a Riemann sum with values of /3 uniformly distributed in 




Figure 3. (Colour online) Fidelity for the QKR (red ) and the pendulum 

(green ) for e = 0.075, ki = 0.6tt and fca = O.Stt. (a) /3 = 0.3216, (b) /3 = 0.2412. 

The data are averaged over 100 kicks in order to get rid of the fast oscillations. 




Figure 4. (Colour online) Comparison between the fidelity for the perturbative 
approach (dashed line) and the QKR (solid line), e = 0.05, fci — O.Gtt, k2 — O.Stt, 
A/3 = 0.06, /3i = 0.06 (black, right), = 0.14 (red, middle), = 0.22 (green, left). 



[Pi, Pi + A/3]. It has been checked that A^^ should be of order of a few thousand to 
get a reasonably good approximation for the integral in (j9]). 

For the boundary of the phase space cell the same criterion as used for single rotors 
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e 


range 


for ensemble 






theoretical upper bound 


0.1 


0.10 - 


-0.16 < /3i,/32 


< 0.36 


-0.37 


0.34 


0.05 


0.12 - 


-0.16 < /3i,/32 


< 0.39 


- 0.41 


0.39 


0.01 


0.12 - 


-0.16 < /3i,/32 


< 0.45 


-0.47 


0.45 



Table 1. Range of validity of the pendulum approximation. For simplicity we defined 
/32 = A/3 + /3i. 




Figure 5. (Colour online) Fidelity of a few ensembles of /3— rotors with e ~ 0.05, 

A/3 = 0.06, ki = 0.67r, fca = O.Stt, /3i = 0.08 (black ), [3i = 0.12 (red ), 

/3i = 0.16 (green ), /3i = 0.2 (blue ) and /3i = 0.24 (magenta —••—). In 

(a) the original data and in (b) the rescaled data are shown. The scaling factors are 
according to figure [6l 



is applied. In Fig. H] a few ensembles are shown. For the measure of correspondence 
we compare the widths of the first few pseudo-oscillations and also demand a deviation 
of less than 10% here. The onset of the island behaviour near to a resonance leads 
to small peaks in the fidelity as described in [S]. Therefore the occurrence of these 
peaks defines the critical value near the resonance island. The intervals for several e are 
summarised in Table [H The upper bound is described quite well by the estimate due 
to the pendulum approximation (|25l) and only fails for the largest e presented, i.e. far 
from the semiclassical regime. 

5. Scaling of ensembles 

In this section we perform additional numerical analysis over a larger range of 
parameters. This will confirm the range of validity of our perturbative approach. 

The ensembles in Fig. [5^ show a similar fidelity as a function of time except for a 
shift along the time axis. This suggests a rescaling of the time. First choose a reference 
value of /3, say /3ref- Then we claim that we can map the fidelity for another f3i on top 
of the reference value only by rescaling the time. This can be more formally written in 
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0.1 



0.2 

Pi 



(a) 



0.3 0.4 



10 



0.1 




0.1 



0.2 

Pi 



0.3 0.4 



(b) 



Figure 6. (Colour online) Scaling factors for several parameters. The scaling is 
done for two references /3rof — 0.16 (upper curves) and /3ref = 0.24 (lower curves), 
fca = O.Stt, = 0.03 (black circles), A/3 = 0.06 (red squares), A^ ^ 0.09 (green 
diamonds), A/3 = 0.12 (blue triangles), [(a)] e = 0.005 Afc ^ 0.2n, and [(b)] e = 0.05 
Afc = O.Itt. 



the following way (/3ref is chosen a priori): 

F(fci, k2, A/3, t) ^ F{ki, k2, (3re{, A/3, a(/3i)-^ t) . (26) 

The fidelity is scaled such that a box surrounding the mean decay can be kept as narrow 
as possible. An example where the curves in Fig. [5^ have been rescaled is presented 
in Fig. |5b. One may notice that the formula (l26ll comes from heuristic and numerical 
observations. It is believed to work as long as the pendulum approximation does for the 
kicked rotor. The latter approximation cannot be controlled in a simple way [9]. That 
is the reason why we will only give here numerical bounds for the validity of ( 126]) . 

This scahng procedure was done for several /3i, A/3, e and ki, whereas ^2 was kept 
the same. The procedure was repeated for another reference (3rei- We expect the curves 
corresponding to different reference values to be parallel, i.e. a different choice of Pj-ei 
just leads to a constant offset. The resulting plots are presented in Fig. El We can see 
that the scaling factor has no strong dependence on any of the parameters shown. 

The same scaling procedure also has been performed for the perturbative data 
computed from flMl) . In Fig. [7] we show the comparison between scaling factors of the 
perturbative result and the QKR for one set of parameters. We observe a reasonable 
agreement. Below /3i ~ 0.12 there is a systematical deviation between Eq. ( 12^ and the 
QKR data. This corresponds to the border of the phase space cell where we already 
expect the pendulum approach to fail. One may conclude that the rescaling procedure 
is valid as long as the pendulum approximation holds. More precisely in our case the 
range of validity is given by the numerical bounds stated above in Sect. 4.2. 
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Pi 

Figure 7. (Colour online) Scaling factors for the perturbative result (black circles) 
and QKR data (red squares), for e = 0.05, ki = 0.67r, fc2 = O.Stt, (i^^i — 0.16 and 
A/3 = 0.06. The dashed line is the border of correspondence between the QKR and 
the perturbative result given in Table [H 

6. Summary and outlook 

We have for the first time studied analytically the quantum fidelity of initial conditions 
corresponding to rotational orbits in the underlying pseudo-classical model. Using 
the pendulum approximation for these orbits, we arrive at our main analytical result 
summarised in Eq. ( 12^ . Although we use a formally somewhat inconsistent expansion 
in the perturbation parameter e, we see that our approximation is rather good when 
including higher orders in the dynamical phases, even if just the lowest order in the 
amplitude of the wave functions is considered. We give clear ranges of applicability 
of our approximation which were tested against numerical simulations of the original 
quantum kicked rotor system. Within these ranges a scaling hypothesis for the temporal 
decay of the fidelity is found which is fulfilled by the original model as well as our 
perturbative results. It would be interesting to set this scaling hypothesis onto firm 
grounds by deriving it from first principles for the here investigated rotational pseudo- 
classical orbits. This task is left for future investigations. 
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We are interested in ttie quantum energy levels of the pendulum Hamiltonian (JTj). It 
is worth reminding that k = k\e\ where |e| is our effective Planck constant. We are 
interested in the regime of small |e|. Our procedure is the following: assume first that 
|e| is a fixed small quantity. Then write perturbative expansions in k, which are valid 
whenever k = 0{\e\). At the very end we will write more explicitly k = k\e\ to derive 
e— semiclassical results. 

The classical action for the Hamiltonian ([7]) along a trajectory from 6i to at a 
fixed energy E is: 



2{E -kcos(fi)d(p - ^{9 - Oi 



(A.l) 



where ^ = sgn(e)(— vr + r/3) = .^o + The energy level can be well described using a 
WKB-like approach in the regime |e| — t- by: 



2n 



2{E^ — k COS (p)d(p — 2^Tc = 2Tcm\e\ . 



(A.2) 



Here the Maslov index is as the particle lives on a ring, hence never meets any 
boundary. The quantization condition f lA.2p with integers m can be rewritten in a 
more efficient way as 



4j2{E^-k)E\u 



2k 



Et 



27r [(m + /3)|e|+eo] 



where E(k) is the Legendre complete elliptic integral: 

'•7r/2 



Jo 



k2 sin(t)2dt . 



(A.3) 



(A.4) 



Eq. ( lA.SIl is the starting point of our perturbation expansion. We assume from now 
that ^ 7^ 0. For the left hand side of flA.3|l the Taylor expansion of is used, see e.g. 
8.114 (1) p.853 in [29]: 



^ ' 2 



K 



64 256 16384 



(A.5) 



One gets: 




+ 



9fc3 241A;^ 



2Ei IGE^J 32E^3 1024^^4 
The eigenenergy E'^ is assumed to have the following form: 



E!i^ 



Et 



so 



(A.6) 



(A.7) 
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Identifying both parts of (lA.Sp using (1A.6P leads to the following results: 




(A.8) 



(m + 



(A.9) 



2^ 



(A.IO) 



4 eo 64 eo 



Inserting f[0|) . (TOj) and ([OOj) into (TOl) gives (fTTj). 
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